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The parameters obtained from the theoretical analysis of the single photon spectra observed by 
the WA98 collaboration at SPS energies have been used to evaluate the two photon correlation 
functions. The single photon spectra and the two photon correlations at RHIC energies have also 
been evaluated, taking into account the effects of the possible spectral change of hadrons in a thermal 
bath. We find that the ratio Rside/Rout ~ 1 for SPS and Rside/Rout < 1 for RHIC energy. 
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I. INTRODUCTION 



Two-particle intensity interferometry along with 
the analysis of single-particle spectra has been 
widely used for a quantitative characterization of 
the hot zone created by the collision of two heavy 
ions at high energies. The method was first ap- 
plied to measure the angular diameter of stars and 
other astronomical objects using the measurements 
based on two-photon correlations known as Hanbury 
Brown-Twiss (HBT) correlations 0. Later in the 
field of particle physics, two-pion correlations were 
used to obtain the spatial extent of the fireball in 
proton - anti-proton reactions 2]. In high energy 
heavy ion collisions two-particle {e.g. pions, kaons 
etc) correlations have been extensively studied both 
experimentally H, ^ |^ and theoretically [E 13 
obtain direct information about the size, shape and 
dynamics of the source at freeze-out. This is usu- 
ally done via selection of transverse momentum and 
rapidity of the correlated particles. Compared to 
two-particle correlation, the three-particle correla- 
tion has been shown to provide additional informa- 
tion on chaoticity and asymmetry of the source emis- 
sion 0. Recently, it has been argued that spatial 
separation between sources of different species can 
be extracted from non-identical particle correlation 
functions, hence these can provide an independent 
cross-check of the transverse flow prescription 0. 
One of the major limitations of carrying out the 
correlation studies with hadrons appearing at the fi- 
nal state is that, the information about the possible 
early dense state of matter is diluted or lost through 
re-scattering. Although, some calculations suggest 
that the ratio of the radius of the source along the 
direction of the total transverse momentum (Kt) 
of the two detected particles to the source radius 
perpendicular to Kt and the beam direction, as a 
function of transverse momentum can provide useful 
information 10]. However, the results at relativistic 
heavy ion collider (RHIC) where quark-gluon plasma 



(QGP) is expected to be formed does not show any 
encouraging result on the above signal 0, . Inves- 
tigations to understand the partonic effects on the 
interferometry with final state particles (pions) at 
RHIC 01 are being pursued. HBT interferometry 
as a sensitive probe of the QCD equation of state 
and hence formation of QGP has been discussed 
in [Hill m. It has been argued [H H that in 
contrast to hadrons, two-particle intensity interfer- 
ometry of photons which are produced throughout 
the space-time evolution of the reaction and which 
suffer almost no interactions with the surrounding 
medium can provide information on the the history 
of the evolution of the hot matter created in heavy 
ion collisions. 

From the experimental point of view, photon 
interferometry encounters considerable difficulties 
compared to hadron interferometry due to small 
yield of direct photons from the early hot and 
dense region of the matter and the associated large 
background of photons primarily from the electro- 
magnetic decay processes of hadrons at freeze-out 
[Tsl l. However, recent calculations demonstrate that 
it is still possible to experimentally filter out the 
correlations of photons from the hot and dense zone 
from those arising due to the residual correlations of 
decay photons. This can be done by studying pho- 
ton interferometry as a function of invariant relative 
momentum of the two photons [l^ . One expects the 
contribution from the photons coming from the early 
stage of the reactions to contribute predominantly to 
the region of small invariant relative momentum. 

The aim of the present work is to first analyze the 
single photon spectra observed by the WA98 collab- 
oration [isj in Pb -I- Pb interactions at SPS ener- 
gies. In the next step we evaluate the two photon 
correlations with the initial conditions, equation of 
states (EOS) and freeze-out conditions which repro- 
duce the photon spectra mentioned above. The EOS 
used here ^3 contains all hadrons upto mass ^ 2.5 
GeV The freeze-out conditions are fixed from 
the study of hadronic spectra measured by NA49 col- 
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laboration 19]. The effects of the possible spectral 
modifications of hadrons in a thermal bath on the 
single photon spectra and two photons correlations 
have been included. The results at RHIC energy 
for the single photon distributions and two photon 
correlations arc also presented. 

It is necessary to point out here that our calcu- 
lation contains some significant improvements over 
the previous works 13, 14]. The correlation func- 
tion evaluated here (i) contain contributions from 
two loop calculations of photon emission rate [20|. 
(ii) emission rate from hadronic matter contain in- 
medium effects, (iii) the EOS contains all hadronic 
degrees of freedom upto mass ~ 2.5 GeV and (iv) 
the initial conditions, freeze-out conditions and the 
EOS reproduces the WA98 single photon spectra |Q 
and NA49 hadronic spectra. Moreover, in Ref. [Tj] 
where the correct definition of the BEG was used, 
the space time evolution was only one-dimensional. 

The paper is organized as follows. In the next 
section we give a general discussion on the correla- 
tion function and the associated kinematics. This is 
followed by the section which deals with the space 
time evolution. Section IV, deals with the results. 
First we present the results for SPS energies. Then 



we discuss the single photon spectra and two pho- 
ton correlations at RHIC energies. In section V we 
summarize our findings. 



II. CORRELATION FUNCTION 

The Bose-Einstein correlation (BEG) function for 
two identical particles is defined as. 



P2{kuk2) 

Pi{h)Pi{k2) 
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where ki is the three momentum of the particle i and 
Pi{ki) and P2{ki,k2) represent the one- and two- 
particle inclusive photon spectra respectively. These 
are defined as, 



P^{k) = / d'^x uj{x,k) 



(2) 



and 
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P2{ki,k2) = Pi{ki)Pi{k2) + / d'xi<fx2 oj{xi,K) uj{x2,K) cos(Ax''AA:^) 
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where, K = {ki + fc2)/2, Afc^ = fci^ - = x 
and k are the four-coordinates of position and mo- 
mentum respectively and Lo{x,k) is the source func- 
tion, which defines the average number of particles 
with four-momentum k emitted from a source ele- 



ment centered at the space-time point x. uj(x, k) is 
actually the thermal emission rate of photons per 
unit four volume. For the QGP phase at a temper- 
ature T this is given by (23, 12^ > 
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where, a — 1/137, is the strong coupling con- 
stant, Jt = 4.45 and Jl = -4.26. Note that C2 is 
independent of the strong coupling constant. For 
photon emission from a hot hadronic gas, we have 
considered a host of reaction processes involving tt, 
p, 77 and uj mesons as well as from the decay of the p 
and u! ^] . For the details on the evaluations of the 
photon emission rate form a hadronic gas we refer 
to [23, 13] . The effects of the intermediary ai meson 



have also been considered in evaluating the photon 
emission rate HM Hfil . 



We shall be presenting the results in terms of dif- 
ference in rapidity (Ay) of the two photons, outward 
(lout), side- ward (qside) and invariant momentum 
differences (gini,) of the two photons and these are 
defined as. 
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where, = kiT — k2T, Kt = (/ciT + fc2T)/2 with the 
subscript T indicating the transverse component, yi 
is the rapidity and ipi's are the angles made by 
with the X-axis. One also defines another kinematic 
variable qiong as 



qiong = kiz - k2z 

= kiT sinh(2/i) - k2T sinh(2/2) 



(8) 



It may be mentioned that the BEG function has 
values 1 < C2(fci,fc2) < 2 for a chaotic source. 
These bounds are from quantum statistics. While 
the radius corresponding to qside {Rside) is closely 
related to the transverse size of the system, the ra- 
dius corresponding to qout {Rout) measures both the 
transverse size and duration of particle emission 6, 
01 ■ So studying Rside/ Rout wiU indicate about the 
duration of particle emission [27i 12^ l29l |. These 
source dimensions can be obtained by parametrizing 
the calculated correlation function with the empiri- 
cal Gaussian form 

^2 = 1+ ^^Vi^ R^outqlimt " Rsideqlide ~ Rlongqiong) 

f9) 

A gross idea of the source size can also be inferred 
from Ri„„ which is defined as 



C2 = 1 + eKY>{-Rlnv% 
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III. SPACE-TIME EVOLUTION 



To compare transverse momentum spectra of pho- 
ton with the experiments we have to convolute the 
static emission rate (fixed temperature, e.g. given 
by eq. ^ by the space-time evolution from the for- 
mation to the freeze-out state. This is done solving 
the relativistic hydrodynamical equations in (3-1-1) 



dimension [30j with boost invariance along the longi- 
tudinal direction . One then needs to specify the 
following: the initial temperature and time when the 
hot system reaches a state of thermodynamic equi- 
librium, the EOS which guides the rate of expan- 
sion/cooling, and the freeze-out temperature when 
the system decouples into free-streaming hadrons. 
We will not discuss here the details on the EOS, ini- 
tial conditions and the modification of the hadronic 
masses in a thermal bath because these topics had 
been discussed in our earlier works in great detail, 
which will be appropriately referred below. For SPS 
conditions we will assume a hot hadronic gas in the 
initial state with temperature dependent masses of 
the constituent hadrons which expands and cools till 
freeze-out. The two-photon correlations at SPS en- 
ergies will also be given for a possible phase transi- 
tion scenario i.e. when QGP is formed at the initial 
state. The EOS used here is given in The 
initial conditions are similar to ref. [s^ Isjj- The 
variation of hadronic masses (except pseudoscalar) 
with temperature follows universal scaling law pro- 
posed by Brown and Rho (BR) j^l- It is necessary 
to point out that a substantial amount of literature 
is devoted to the calculation of effective masses of 
hadrons in the medium. Our principal motivation 
for choosing the BR conjecture is that this has been 
used successfully to explain the WA98 photon spec- 
tra Hi and GERES/NA45 [13 low mass dilcpton 
spectra 1^3 • For higher collision energies, i.e. 
at RHIG we will consider a situation where QGP is 
formed in the initial state, then the quark matter 
evolves with time to a hadronic phase via an inter- 
mediate mixed phase in a first order phase transi- 
tion scenario. The mixed phase is a mixture of both 
quark matter and hadronic matter, with the fraction 
of quark matter decreasing with time to zero when 
the phase transition is complete. The hot hadronic 
gas then expands till the system freezes-out. The 
initial condition for the SPS and RHIG energies in 
terms of the initial temperature (T^) is set from the 
number of particles per unit rapidity at the mid- 
rapidity region for those energies according to the 
following equation. 
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o-k — 7r^(7fc/90 is determined by the statistical degen- 
eracy {gk) of the system formed after the collision. 
Taking the particle multiplicity per unit rapidity at 
mid rapidity to be 700 for SPS, we get T, = 200 
MeV with the formation time (r^) taken as 1.0 fm/c. 
With a multiplicity of 1100 for Au-|-Au collisions at 
RHIG, we get Ti = 264 MeV for an initial time of 
0.6 fm/c (33. The critical temperature (Tc) is taken 
to be 170 MeV [l^ here. The freeze-out tempera- 
ture (T/) is taken to be 120 MeV for both SPS and 
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RHIC energies; a value that describes the transverse 
momentum distributions of hadrons produced |40l | . 
Let us now turn to the EOS which plays a central 
role in the space time evolution, we have consid- 
ered. For the hadronic phase the EOS corresponds 
to that for the hadronic gas with particles of mass 
up to 2.5 GeV and includes the effects of non-zero 
widths of various mesonic and hadronic degrees of 
freedom 0. The velocity of sound (Cg) correspond- 
ing to this EOS at freeze-out is about 0.18 [13 (the 
value corresponding to a free gas of massless par- 
ticles being 0.33). We emphasize that the medium 
modifications of the constituent hadrons play a non 
trivial role in the estimation of the quantities men- 
tioned above. Of special mention is the velocity of 
sound and the estimation of the initial temperature 
for the hot hadronic gas considered at SPS. These 
have been exhaustively dealt with in a number of 
our earlier works 0, Is^, E3 and we do not re- 
peat them here. With these inputs we have per- 
formed a (3-1-1 )-dimensional hydrodynamic expan- 
sion. The solution of the hydrodynamic equations 
has been used to evaluate the space time integration 
involved in Eqs. [21and|21 



IV. RESULTS 

In this section we present the results of two- 
photon interferometry. We shall first present the re- 
sults for SPS energies and then give our predictions 
for RHIC. In both the cases, we will first start with 
the single photon spectra and then present the two- 
photon correlation as a function of Ay, qout, Iside 
and Qinv At the SPS energies we present the two- 
photon correlation results using the inputs which re- 
produces the the measured single photon spectra |0| 
as mentioned before. For RHIC our predictions are 
based on the same model as that explains the SPS 
data, but for different initial conditions as expected 
at RHIC energies. 

Fig-ffl shows the single photon yield measured by 
the WA98 experiment at CERN SPS, as a function 
of transverse momentum. The contributions from 
hard scattering of partons (hard photons) and from a 
hot hadronic gas (thermal photons) have been shown 
separately. The hard photon contribution has been 
normalized to reproduce the scaled p-|-C data 
as well as the scaled p-p data [ia (p-p data is not 
shown here. Please see ref. fg^). The experimen- 
tal results of WA98 collaboration show some "ex- 
cess" photons over the hard photons for 1.5 < pr 
(GeV) < 2.5 GeV, which we argue to originate from 
a thermal source either of QGP or hot hadronic gas 
of initial temperature 200 MeV (see also [l^, for 
details). In the present work we have evaluated the 
thermal photon from a hot hadronic gas with their 
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FIG. 1: Single photon spectra at SPS 
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FIG. 2: Correlation function C2 as a function of Ay, 
Qout, Qside and qinv for Pb+Pb collisions at 158 AGeV at 
SPS 



mass reduced according to the universal scaling sce- 
nario |34| . The combined results of hard and thermal 
photons is found to explain the data reasonably well. 

From the expressions given in Eqs. |31|SlE|one can 
see that several combinations can be made in the 
variables y, ip and kx to present the results for two- 
photon correlation. For simplicity, we will present 
all two-photon correlation results for direct photons 
with momentum around 2 GeV/c since it lies in 
the range of momentum where excess direct pho- 
tons have been observed at SPS energies. We take 
■02 = and y2 = for all cases. We vary and yi 
wherever necessary. 

Figs. 121 a-d) show the variation of the correlation 
strength (C2) as a function of Ay, qout, Qside and qinv 
for various phases (sum = QGP-|-Mixed-|-hadrons). 
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FIG. 3: Correlation function, C2 as a function of Ay, 
Qout, Qside and qinv for Pb + Pb collisions at 158 AGeV 
at SPS. Solid (dashed) line indicate results for QGP 
(Mixed) phase and dot-dashed (dotted) line represent 
correlation function for hadronic (sum) phase. 



The HBT radii of the evolving hot hadronic mat- 
ter (denoted by Hadron* indicating that the masses 
of the hadrons are modified in the medium accord- 
ing to BR scaling) is shown in Table 1. Please note 
that to get the quantity Riong in the longitudinally 
co-moving system (LCMS)of reference one should 
multiply the numbers given in Table I by cosh^yn) 
where yK is the rapidity corresponding to the mo- 
mentum K defined in section II. We would like to 
mention here that the HBT radii give the length of 
homogeneity of the source and this is equal to the 
geometric size if the source is static. However, for a 
dynamic source e.g. the system formed after ultra- 
relativistic heavy ion collisions, the HBT radii is 
smaller than the geometric size (see |i28».i43t.i44i45j] ). 

Our earlier analysis js^ shows that the WA98 sin- 
gle photon spectra can be explained by two different 
kind of initial states: (i) hot hadronic gas, with the 
masses of the hadrons vary according to BR scal- 
ing hypothesis or (ii) QGP is created initially. It 
is rather difficult to distinguish between the two at 
present. Therefore, in the following we show the 
results on two-photon correlation functions for case 
(ii) also. In Fig. O the variation of the two-photon 
correlations is depicted in a scenario when QGP is 
formed initially at SPS. The initial conditions here 
are similar to the that of Ref. HI. The HBT di- 
mensions extracted from these correlation functions 
are shown in Table I. The difference in the corre- 
lation function as a function of qout for the QGP 
and the mixed phase is negligible. When plotted 



as a function of qside (Fig- Eb) the correlation func- 
tion for the QGP, mixed phase and the sum almost 
overlap. The width of the correlation function for 
the hadronic phase is the largest as compared to the 
other phases. The HBT dimensions satisfy the re- 
lation Rside/Rout ^ 1 for the correlation functions 
denoted by 'sum' in the Table I. This is irrespec- 
tive of the formation of the QGP or hot hadronic 
gas with the mass of the hadrons (except pseudo- 
scalar) approaching zero at the initial temperature. 
Although the later scenario is closely related to the 
chiral/dcconfinement phase transition it is very dif- 
ficult at this stage to draw a firm conclusion. 



Fig. ^ shows our predictions for the single photon 
spectra at RHIC energies. Here the thermal pho- 
tons contain contribution from QGP, mixed and the 
hadronic phases. The results show a clear domi- 
nance of the thermal photon over the hard photons 
for < 3 GeV. We calculate the two-photon cor- 
relation at RHIC energies with the inputs that are 
used to evaluate the single photon spectra as a func- 
tion of Ay, Qout, Qside and qtnv 
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FIG. 5: Correlation function C2 as a function of Ay, 
Qout, Qside and Qinv ioT Au+Au collisious at 200 AGeV 
at RHIC. Solid (dashed) line indicate results for QGP 
(Mixed) phase and dot-dashed (dotted) line represent 
correlation function for hadronic (sum) phase. 







Rinv 


Rout 


Rside 


Rlong 


SPS 


Hadron* 


3.7 


3.6 


3.7 


0.5 




QGP 


3.5 


3.4 


3.5 


0.37 


SPS 


Mixed 


3.5 


3.4 


3.5 


0.8 




Hadron 


3.0 


4.5 


3.0 


1.5 




Sum 


3.5 


3.6 


3.5 


0.9 




QGP 


3.5 


3.5 


3.5 


0.55 


RHIC 


Mixed 


3 


5.4 


3 


2.0 




Hadron 


2.7 


6.8 


2.7 


2.3 




Sum 


3 


5.5 


3 


1.6 



Table 1 : Values of the various parameters (in fm) of 
the correlation functions in the forms given in Eqs. 
1^ and [TUl for SPS and RHIC energies at a transverse 
momentum 2 GeV. 



Fig. [S^ shows the variation of C2 as a function of 
Ay, for two photons with momentum of 2 GeV/c. 
We have taken ipi — ip2 — for simplicity. The 
contribution from all the three phases are shown. 
We observe that the width of the correlation func- 
tion is largest for the QGP phase followed by that 
from the mixed phase and then the hadronic phase. 
The inverse will reflect the corresponding HBT radii 
of the source. Fig.[SjD shows the variation of C2 with 
Qout, for yi = 2/2 = 0, V^i = -02 =0 and kix 
— 2 GeV. The contribution of all three phases are 
shown. One observes that the C2 variation is sim- 
ilar for the three phases and is similar to that for 
the case with Ay but the widths are considerably 
smaller. Fig. shows how C2 varies with qside for 
all the three phases as well as the sum. For this 
case we have taken, yi = 2/2 = 0, ?/'2 = and 
kiT — k2T — 2 GeV. The trend in the varia- 
tion of the width for the three phase has reversed. 
The width of C2 distribution for QGP is smallest, 
followed by that for the mixed phase and then the 
hadronic phase. Fig. [SJi shows how C2 varies with 
Qinv for all the three phases as well as the sum. For 
this case we have taken, yi = j/2 = 0, '02 =0 and 
kiT = k2T = 2 GeV. The trend in the variation 
of the width for the three phases is similar to that 
for the case of qside- 



V. SUMMARY 

The two photon correlation functions has been 
evaluated for SPS and RHIC energies. Constraints 
from the experimentally observed single photon 
spectra has been used to evaluate the correla- 
tion functions at SPS energies. Predictions for 
both the single photon spectra and the two pho- 
ton correlation functions have also been given at 
RHIC energies. The values of HBT radii extracted 
from the two-photon correlation functions show that 
Rside/Rout ~ 1 irrespective of the formation of QGP 
or hadronic gas with reduced mass initially for SPS 
and Rside/Rout < 1 for RHIC energy. 
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